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I Abstract. We compute spin-dependent matrix elements for decays of 5-wave quarkonia in lattice 

. NRQCD. They appear to be in approximate agreement with the velocity-scaling rules of NRQCD. 

Q ! PACS: 12.38Ge,13.25Gv 

Q: 

(M ; The effective field tlieory Nonrelativistic QCD (NRQCD) H S H is an elegant 

\ formalism within which to describe heavy-quarkonium physics. It is believed that, 

^ ■ within NRQCD, inclusive cross sections for the production of a quarkonium at large 

\0 '. transverse momentum pj can be written in a factorized form [3]: 

O; a{H)=Y^an{mWn{m). (1) 

^ . 

^ \ where the a„ are perturbatively calculable short-distance coefficients, and the 

■ (0|^^(A)|0) are matrix elements in the vacuum state of four-fermion operators of 

! the form = X^Kn^fiXx \H -\-X){H -\-X\) xif^K^X- Here, H is the quarkonium state, 

I// is the Pauli spinor field that annihilates a heavy quark, x is the Pauli spinor field that 
creates a heavy antiquark, and A is the ultraviolet cutoff of NRQCD. K contains Pauli 
^ . matrices, color matrices, and the covariant derivatives. NRQCD predicts the leading 

scaling behavior of the matrix elements with v, the heavy-quark velocity in the quarko- 
nium rest frame As a consequence of these v-scaling rules, the sum over operator 
^ matrix elements can be regarded as an expansion in powers of v, where ~ 0.3 for 

5^ I charmonium and ^0.1 for bottomonium. A similar factorization formula applies 

to inclusive quarkonium decays [jj], except that the matrix elements are now between 
quarkonium states, rather than vacuum states, and the four-fermion operators have the 

form ffn = W^KnXX^KW- 

At large pr at the Tevatron, the dominant mechanism for 7/i// production is gluon 
fragmentation into a QQ pair, which then evolves nonperturbatively into the J /^f. The 
nonperturbative evolution is described by the NRQCD matrix elements. The NRQCD 
v-scaling rules predict that the spin-flip matrix elements are suppressed by in com- 
parison with the non-spin-flip matrix elements. Therefore, in existing calculations, the 
spin-flip effects are neglected, and the 7/ y/^ is assumed to take on the polarization of the 
fragmenting gluon. Consequently, the J / Xj/ is expected to have a significant transverse 
polarization at large pr (Ref. [^). Surprisingly, the CDF data for the J/ Xjr polarization 
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Js*] show decreasing transverse polarization with increasing pr and disagree with the 
NRQCD prediction [6] in the largest pj bin. 

One can question whether the neglect of spin-flip processes is justified in the case of 
the J /^f- The v-scaling rules predict the leading power of v in a matrix element, but not 
its coefficient. It could happen that the spin-flip matrix elements, although suppressed 
by relative to the non-spin-flip matrix elements, are not actually smaller numerically. 
It has also been suggested that the v-scaling rules themselves may need to be modified 
for charmonium iSESIIUnil. 

Ideally, we would settle these issues by computing the 7/i// production matrix ele- 
ments in a lattice QCD simulation. Unfortunately, it is not known how to formulate the 
computation of production matrix elements in lattice simulations. However, we can in- 
stead use lattice calculations of decay matrix elements to test the validity of using the 
v-scaling rules to estimate the sizes of matrix elements. 

We take for our lattice action the discretized version of the NRQCD action through 
next-to-leading order in v that is given in Ref. We also include improvements 
of relative-order a to the terms in the lattice action that are of leading order in v^, 
and we implement tadpole improvement [13] as described in Ref. II12I1 . Our initial 
computations were carried out on quenched gauge-field configurations, which should 
reproduce the qualitative features of QCD. Our preliminary results are based on 400 
gauge-field configurations on 12^ x 24 lattices at /3 = 5.7. We use the values anic = 0.8 
and amh = 3.15 for the heavy-quark masses in lattice units. We take the parameter n that 
appears in the action of Ref. [12] to be 1 for bottomonium and 4 for charmonium. 

Since, at /3 = 5.7, a = 0.81 GeV^i for charmonium JH and a = 0.73 GeV^i for 
bottomonium lITill . the quarkonium is well contained in the lattice volume. However, 
since the quarkonium radius is approximately l/(mv), which is about 1.2 GeV^^ for 
charmonium and about 0.6 GeV^ bottomonium, the lattice spacing is fairly coarse and 
warrants the implementation of order-a improvements. 

Our preliminary results for 5- wave quarkonium states are shown in Tabled As can be 
seen, the hierarchy of matrix elements that is predicted by v-scaling is observed. How- 
ever, the sizes of the matrix elements decrease faster with increasing powers of v than 
one would expect from the powers of v alone. The NRQCD heavy-quark spin symmetry 
^ predicts that the singlet-to-triplet transition and the triplet-to-singlet transition should 
be in a ratio of 3: 1, up to corrections of order v^. The lattice results agree with this pre- 
diction. Furthermore, the ratios of bottomonium matrix elements to charmonium matrix 
elements for the larger, well-measured matrix elements agree within about a factor of 
two with expectations from v scaling. Taken together, these results suggest that v scaling 
is fairly well obeyed, but that additional factors, beyond the powers v, must be taken into 
account in order to estimate the sizes of the matrix elements accurately. 

As can be seen from Tabled the triplet-to-singlet transition rate is comparable to the 
triplet-to-triplet transition rate. This suggests that, at large pj at the Tevatron, where 
the color-octet, spin-triplet process dominates 5-wave quarkonium production, the r\c 
production rate may be comparable to the y/i// production rate. 

The spin-triplet up-to-longitudinal transition rate is small compared with the spin- 
triplet up-to-up transition rate. That is, the spin-flip matrix elements are suppressed, as 
is expected from v-scaling. Hence, our preliminary results support the prediction of large 
transverse polarization of 7/i/a's produced at large pj at the Tevatron. 



TABLE 1. Quarkonium decay matrix elements. A spin transition labeled "5/ S/' refers to the matrix 
element of the color-octet 5-wave operator of spin Sf in an 5-wave color-singlet state of spin 5,-. The 
matrix elements are normalized to the matrix element of the ^Si color-singlet operator in a ^Si color- 
singlet state. We average over unspecified spins in the quarkonium state and sum over unspecified spins in 
the operators. The column labeled "scaling" gives the v scaling of the normalized matrix element. Here we 
include the color factor 1 / {2Nc) that arises in the free QQ matrix elements, as is suggested in Ref. Iisll . 
The columns labeled "bottom." and "charm." give the lattice results for the normalized matrix elements 
in the lowest-lying S'-wave bottomonium and charmonium states. The columns labeled "bottom, est." and 
"charm, est." give the numerical values of the velocity-scaling estimates in the "scaling" column, taking 
=0.1 for bottomonium and = 0.3 for charmonium. The quoted uncertainties are statistical only. 
Systematic uncertainties could be considerably larger. 
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